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Abstract: We study the Navier-Stokes equations in dimension 3 (NS3D) driven 
by a noise which is white in time. We establish that if the noise is at same time 
sufficiently smooth and non degenerate in space, then the weak solutions converge 
exponentially fast to equilibrium. 

We use a coupling method. The arguments used in dimension two do not apply 
since, as is well known, uniqueness is an open problem for NS3D. New ideas are 
introduced. Note however that many simplifications appears since we work with 
non degenerate noises. 
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Introduction 


We are concerned with the stochastic Navier-Stokes equations on a three di¬ 
mensional bounded domain (NS3D) with Dirichlet boundary conditions. These 
equations describe the time evolution of an incompressible fluid subjected to a 
determinist and a random exterior force and are given by 


( 0 . 1 ) 


dX + v(- A)A dt + (A, V)X dt + Xpdt = </>(X)dW + f dt, 
(divl)(t,() = 0, for £ G D, t > 0, 

X(t,£) = 0, for £ £ dD, t> 0, 

*((U) = *o(0, for £ € D. 


Here D is an open bounded domain of R 3 with smooth boundary dD or D = (0, l) 3 . 
We have denoted by X the velocity, by p the pressure and by v the viscosity. The 
external force field acting on the fluid is the sum of a random force field of white 
noise type (f>(X)d,W and a determinist one f dt. 

In the deterministic case {<j> = 0), there exists a global weak solution (in the 
PDE sense) of ED when xq is square integrable, but uniqueness of such solution 
is not known. On another hand, there exists a unique local strong solution when xq 
is smooth, but global existence is an open problem (see m for a survey on these 
questions). 
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In the stochastic case, there exists a global weak solution of the martingale 
problem, but pathwise uniqueness or uniqueness in law remain open problems (see 
da for a survey on the stochastic case). 

The main result of the present article is to establish that, if <f> is at the same 
time sufficiently smooth and non degenerate, then the solutions converge exponen¬ 
tially fast to equilibrium. More precisely, given a solution, there exists a stationary 
solution (which might depends on the given solution), such that the total varia¬ 
tion distance between the laws of the given solution and of the stationary solution 
converges to zero exponentially fast. 

Due to the lack of uniqueness, it is not straightforward to define a Markov 
evolution associated to ED- Some recent progress have been obtained in this 
direction. In El, ID, under conditions on <f> and / very similar to ours, it is shown 
that every solution of ED limit of Galerkin approximations verify the weak Markov 
property. Uniqueness in law is not known but we think that this result is a step 
in this direction. Our result combined with this result implies that the transition 
semi-group constructed in j3| is exponentially mixing. 

Note also that recently, a Markov selection argument has allowed the construc¬ 
tion of a Markov evolution in na- Our result does not directly apply since we only 
consider solutions which are limit of Galerkin approximations. However, suitable 
modifications of our proof might imply that under suitable assumptions on the 
noise, the Markov semi-group constructed in m is also exponentially mixing. 

Our proof relies on coupling arguments. These have been introduced recently in 
the context of stochastic partial differential equations by several authors (see hd, 
ED, El, ED, ED, ED, ED, ED and (52]). The aim was to prove exponential 
mixing for degenerate noise. It was previously observed that the degeneracy of the 
noise on some subspace could be compensated by dissipativity arguments m, 0, 
ED More recently, highly degenerate noise noises have been considered in CD, 

ED 

In all these articles, global well posedness of the stochastic equation is strongly 
used in many places of the proof. As already mentioned, this is not the case for 
the three dimensional Navier-Stokes equations considered here. Thus substantial 
changes in the proof have to be introduced. However, we require that the noise 
is sufficiently non degenerate and many difficulties of the above mentioned articles 
disappear. 

The main idea is that coupling of solutions can be achieved for initial data which 
are small in a sufficiently smooth norm. A coupling satisfying good properties is 
constructed thanks to the Bismut-Elworthy-Li formula. Another important ingre¬ 
dient in our proof is that any weak solution enters a small ball in the smooth norm 
and that the time of entering in this ball admits an exponential moment. We 
overcome the lack of uniqueness of solutions by working with Galerkin approxima¬ 
tions. We prove exponential mixing for these with constants which are controlled 
uniformly. Taking the limit, we obtain our result for solutions which are limit of 
Galerkin approximations. 
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1. Preliminaries and main result 


1.1. Weak solutions. 

Here C(Ki;K 2 ) (resp B 2 {K\\ K 2 )) denotes the space of bounded (resp Hilbert- 
Schmidt) linear operators from the Hilbert space K\ to K 2 . 

We denote by |-| and (•,•) the norm and the inner product of L 2 (Z5;R 3 ) and 
by |-| the norm of L P (D; R 3 ). Recall now the definition of the Sobolev spaces 
H p (D;R 3 ) for p G N 

j HP{D;R 3 ) = {X G L 2 {D-R 3 ) \ d a X G L 2 {D- R 3 ) for [a|<p}, 

1 \ X \hp = T / \ a \< p \ d ‘* X \ 2 ■ 

It is well known that (22 P (Z); R 3 ), is a Hilbert space. The Sobolev space 

Hq{D]R 3 ) is the closure of the space of smooth functions on D with compact 
support by |-| H i. Setting ||A'|| = |VA|, we obtain that ||-|| and are two 

equivalent norms on Hq(D;R 3 ) and that (Hq(D;R 3 ), ||-||) is a Hilbert space. 

Let H and V be the closure of the space of smooth functions on D with compact 
support and free divergence for the norm |-| and ||-||, respectively. 

Let 7r be the orthogonal projection in L 2 (Z5;R 3 ) onto the space H. We set 

A = 7T (—A), Z5(A) =Vn H 2 (D; R 3 ), B(u, v) = n ((u, V)u) and B(u) = B(u, u). 
Let us recall the following useful identities 

{ (B(u,v),v) = 0, u, v G V, 

( B(u,v),w ) = —( B(u,w),v ), u,v,wGV. 

As is classical, we get rid of the pressure and rewrite problem gh> in the form 


{ dX + vAXdt + B{X)dt = <l>(X)dW + f dt, 

X(0) = x 0 , 


where W is a cylindrical Wiener process on H and with a slight abuse of notations, 
we have denoted by the same symbols the projections of 4> and /. 

It is well-known that (A, T>(A)) is a self-adjoint operator with discrete spectrum. 
See ?2Ji El- We consider ( e n ) n an eigenbasis of H associated to the increasing 
sequence (/x ra )„ of eigenvalues of (A,25(A)). It will be convenient to use the frac- 
tionnal power (A S ,2?(A S )) of the operator (A, 25(A)) for s£l 


| T>{A S ) 

j A S X 


E OO o 

n=l ^n X n e n 


W°° Li 
Z^in-1 r’n 


2s 'x n \ 2 < OO 


}• 


where X = 


— i '^n z 'n • 


We set for any s G R 

\\X\\ a = \A%X\, H S =25(A5). 

It is obvious that (H s , ||-|| s ) is a Hilbert space, that (Ho, |HI 0 ) = {H, | j) and that 
(Hi, H'l^) = ( V ,, ||-||). Moreover, recall that, thanks to the regularity theory of the 
Stokes operator, H s is a closed subspace of H S (D, R 3 ) and ||-|| s is equivalent to the 
usual norm of H s (D; R 3 ) when D is an open bounded domain of M 3 with smooth 
boundary dD. When D = (0, l) 3 , it remains true for s < 2. 
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Let us define 

* = L^ c (R+;H)nLl c (R+-V)nC(R+;M s ), 

< W = C(K + ; H_ 2 ), 

£2* = x x w, 

where s is any fixed negative number. Remark that the definition of X is not 
depending on s < 0. Let X* (resp W*) be the projector £2* —► X (resp £2* —> W). 
The space £2* is endowed with its borelian cr-algebra T* and with (X t *) t>0 the 
filtration generated by (X*, W*). 

Recall that W is said to be a (Pt) t -cylindrical Wiener process on H if W is 
(J r t) t -adapted, if W(t + •) — W{t) is independant of T t for any t > 0 and if W is 
a cylindrical Wiener process on H. Let if be a Polish space. We denote by P(E) 
the set of probability measure on E endowed with the borelian er algebra. 

Definition 1.1 (Weak solutions). A probability measure Pa on (£2*,X*) is said to 
be a weak solution of o with initial law A £ P(H ) if the three following properties 
hold. 

i) The law o/X*(0) under Pa is A. 

ii) The process W * is a (Pf)t-cylindrical Wiener process on H under Pa- 

iii) We have Pa -almost surely 

^ ^ (X*(2),'0) + ^/o(X*(s),AV))ds + / o t (S(X*(s)),^)ds 

= (X*(0),V>) +t(f,ip) + / 0 t (V>,<^(X*(s))dW*(s)), 

for any t £ R + and any if smooth mapping on D with compact support and 
divergence zero. 

When the initial value A is not specified, iro is the initial value of the weak solution 
P Xo (i.e. A is equal to S Xo the Dirac mass at point xq). 

These solutions are weak in both probability and PDE sense. On the one hand, 
these are solutions in law. Existence of solutions in law does not imply that, given 
a Wiener process W and an initial condition xq, there exist a solution X associated 
to W and xo- On the other hand, these solutions live in H and it is not known if 
they live in Hi. This latter fact causes many problems when trying to apply Ito 
Formula on X(X*(f)) when F is a smooth mapping. Actually, we do not know if 
we are allowed to apply it. 

That is the reason why we do not consider any weak solution but only those which 
are limit in distribution of solutions of Galerkin approximations of HU- More 
precisely, for any IV £ N, we denote by Pn the eigenprojector of A associated to 
the first N eigenvalues. Let (£2, T, P) be a probability space and W be a cylindrical 
Wiener process on H for P. We consider the following approximation of o 

{ dXpf + izAXjydt + P^B{Xisr)dt = Piy(j)(X]\f)dW + Pjyfdt, 

Xat(0) = Pnxq- 

In order to have existence of a weak solution, we use the following assumption. 

Hypothesis 1.2. The mapping </> is bounded Lipschitz H —> £2 ( H ; Hi) and f £ H. 
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We set 

2 Iff 

B l = SU P \<t>( X )\c 2 {H-Mi) + —■ 

x£H Z'M i 

It is easily shown that, given xo £ H , o has a unique solution Xjf = Xn(-,xo). 
Proceeding as in EH, we can see that the laws (P^)jv of (X;v(-, £o)j W) are tight in 
a well chosen functional space. Then, for a subsequence ( Nk)k , (Xjv fc , W) converges 
in law to P Xo a weak solution of O- Hence we have existence of the weak solutions 
of ltT~Tl) . but uniqueness remains an open problem. 

Remark 1.3. We only consider weak solutions constructed in that way. This 
allows to make some computations and to obtain many estimates. For instance, 
when trying to estimate the L 2 -norm of X*(t) under a weak solution Px 0 , we would 
like to apply the Ito Formula on |X*| . This would give 

d |X*| 2 + 2is\\X4 2 dt = 2(X*,<t>(X*)dW*)+2(f,X*)dt + \4>{X*{t))\ 2 C2{H . H] dt. 

Integrating and taking the expectation, we would deduce that, if f = 0 and (f> con¬ 
stant, 

Exo (V*(X)| 2 + 2 " J Q ll*.(«)H 2 dt) = ko| 2 + t \4>\ 2 c 2{h . h) ■ 

Unfortunately, those computations are not allowed. However, analogous computa¬ 
tions are valid if we replace Px 0 by P^, which yields 

E |jWv(f)| 2 + 2 v J \\X N (s)\\~ dt) = \PnXq\ 2 +t\PN4>\ 2 £ 2 ( H . H ) ■ 

Then, we take the limit and we infer from Fatou Lemma and from the semi¬ 
continuity of |-|, ||-|| in H a that 

Ex 0 (V*(t)| 2 + 2 " J q ll*.(«)ll 2 dt) < |z 0 | 2 + t \<t >\ 2 C2{H . H) , 

provided f = 0 and <f constant and provided P Xo is limit in distribution of solutions 

of 11.311 . 

Let P' and Y be a probability measure and a random variable on (12*, J 7 *), 
respectively. The distribution XV 00 denotes the law of Y under P'. 

A weak solution P M with initial law p is said to be stationary if, for any t > 0, 
p is equal to Vp (X*(t)). 

We define 

(ZfVO (*o) = E (W: N (t, xo))) = e^ o (i>(X*(t)), 

where is the expectation associated to P^. 

It is easily shown that Xn(-,xo) verifies the strong Markov property, which 
obviously implies that ('P t JV ) teR + is a Markov transition semi-group on PnH. 

Ito Formula on |Xjv(-,xo)| gives 

d \X N \ 2 + 2u HAT^vll 2 dt = 2 (X N , </>(X N )dW) + 2{X N , f)dt + \P N cf>{X N )\ 2 dt, 

which yields, by applying arithmetico-geometric inequality and Hypothesis 11.21 

(1.4) d\X N \ 2 + v\\X N \\ 2 dt <2{X N ,(j){X N )dW) + cB 1 dt. 
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Integrating and taking the expectation, we obtain 

(1.5) E (W(f)| 2 ) < e-^ |zo| 2 + —Bi. 

Hence, applying the Krylov-Bogoliubov Criterion (see 0J), we obtain that {V^)t 
admits an invariant measure and that every invariant measure has a moment 
of order two in H. Let Xq be a random variable whose law is and which 
is independent of W, then Xjv = X^(-,Xq) is a stationary solution of 11.311 . 
Integrating El, we obtain 

E\X N {t)\ 2 + pE f ||A' A r(s)|| 2 ds<E|X Ar (0)| 2 + cHit. 

Jo 

Since the law of X/v(s) is pn for any s > 0 and since pn admits a moment of order 
2, it follows 

(1.6) [ \\x\\ 2 p N (dx) < -B 1 . 

Jp n h v 

Moreover the laws (P^)jv °f (-Xjv( - , X^), W) are tight in a well chosen functional 
space. Then, for a subsequence (N' k )k, converges in law to a weak station¬ 
ary solution of ED with initial law p (See El for details). We deduce from 11.(il) 
that 

[ INI 2 IJ-{dx) < -Hi, 

Jh v 

which yields (see El) 

(1.7) (X*(f) G Hi) = 1 for any t> 0. 

We do not know if A* (t) G Hi for all t holds P^-almost surely. This would probably 
imply strong uniqueness /r-almost surely. Remark that it is not known in general 
if /i is an invariant measure because, due to the lack of uniqueness, it is not known 
if ED defines a Markov evolution. We will see below that this is the case under 
suitable assumptions. 

1.2. Exponential convergence to equilibrium. 

In the present article, the covariance operator cf> of the noise is assumed to be 
at the same time sufficiently smooth and non degenerate with bounded derivatives. 
More precisely, we use the following assumption. 

Hypothesis 1.4. There exist e > 0 and a family ((f> n ) n of continuous mappings 
H —> R with continuous derivatives such that 

j 4>(x)dW = Jf l ^L 1 4> n (x)e n dW n where W = W n e n , 

\ K 0 = E,T=1 SU P xeH I^OzOI 2 hn +E < OO. 

Moreover there exists K\ such that for any x, rj G H 2 

OO 

■ h\ 2 hi <K1\\V\\1- 

n =1 

For any x G H and N G N, we have <p n (x) > 0 and 

(1.8) K 2 = sup < OO, 
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where 

OO OO 

• h = ^2 < l ) n(x)~ 1 h n e n for h = h n e n . 
n =1 n=0 

For instance, <j> = A~i fulfills Hypothesis ^2] provided s £ (|,3]. 

We set 

Bo = no + ki + K2 + \ f\ 2 ■ 

Remark 1.5 (Additive noise). If the noise is additive, Hypothesis O simplifies. 
Indeed in this case, we do not need to assume that </> and A commute. This requires 
a different but simpler proof of Lem.m.n fTA below. 

Remark 1.6 (Large viscosity). Another situation where we can get rid of the 
assumption that the noise is diagonal is when the viscosity v is sufficiently large. 
The proof is simpler in that case. 

Remark 1.7. It is easily shown that Hypothesis E3 and f £ H imply Hypothesis 
11.21 Therefore, solutions of O are well-defined and, for a subsequence, they 
converge to weak solution of CD- 

The aim of the present article is to establish that, under Hypothesis o and 
under a condition of smallness of ||/|| e , the law of X*(f) under a weak solution P Xo 
converges exponentially fast to equilibrium provided P Xo is limit in distribution of 
solutions of d. 

Before stating our main result, let us recall some definitions. Let E be Polish 
space. The set of all probability measures on E is denoted by V{E). The set of all 
bounded measurable (resp uniformly continuous) maps from E to R is denoted by 
(resp UCb(E; R)). The total variation of a finite real measure A 

on E is given by 

l|A|L r =sup{|A(T)| | T G B(E)} , 

where we denote by B(E) the set of the Borelian subsets of E. 

The main result of the present article is the following. Its proof is given in section 
4 after several preliminary results. 

Theorem 1.8. Assume that Hypothesis na holds. There exists S°, C and 7 > 0 
only depending on <p, D, £ and v such that, for any weak solution Pa with initial 
law A £ V(H) which is limit of solutions of o> there exists a weak stationary 
solution P^ with initial law /i such that 

(1.9) ||9p x (*.(*)) - MlU < Ce (l + J h M 2 A (dx)^j , 

provided ||/|| 2 < <5° and where ||-|| is the total variation norm associated to the 
space H s for s < 0. 

Moreover, for a given Pa, p is unique and P M is limit of solutions of O- 

It is well known that ||-|| is the dual norm of which means that for any 
finite measure A' on H s for s < 0 

l! A 1Lr= SU P / 9(x)\'(dx) , 
lsloo<l -'His 
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where the supremum is taken over g £ UC b (M s ) which verifies |g| < 1. Hence 

G3J is equivalent to 


( 1 . 10 ) 


Ea (g(X*(t))) - / g(x)p(dx) 
J H 


<C\g\ c 


X (dx) 


/ H 


for any g £ UC b (M. s ). 


Remark 1.9 (Topology associated to the total variation norm). Remark that if A' 
is a finite measure of H So , then the value of the total variation norm of X 1 associated 
to the space H s is not depending of the value of s < so- 

Hence, since Vp x (X*(t)) is a probability measure on H then ( 11.011 (resp ( 11 . 101 ) ) 
remains true when ||-|| is the total variation norm associated to the space H (resp 
for any g £ B b (H;R)). 

Moreover, we see below that, under suitable assumptions, if X is a probability 
measure on H 2; then T>p x (X*(t)) is still a probability measure on H 2 . It follows 
that 03 (resp 11.101 ) remains true when ||-|| is associated to H 2 (resp for any 
g £ H{,(H 2 ; K)j. 


Our method is not influenced by the size of the viscosity v. Then, for simplicity 
in the redaction, we now assume that v = 1. 


1.3. Markov evolution. 

Here, we take into account the results of 0 , m and we rewrite Theorem Ol 
This section is not necessary in the understanding of the proof of Theorem II. 81 
Let (N(,)k be an increasing sequence of integer. In 0 , 0 , it is established that 
it is possible to extract a subsequence ( Nk)k of ( N' k )k such that, for any xq £ H 2 , 
converges in distribution to a weak solution IP Xo of <0> provided the following 
assumption holds. 

Hypothesis 1.10. There exist e,S > 0 such that the mapping cj) is bounded in 
C 2 (if;IHIi +e ). Moreover, for any x, ker <f(x) = {0} and there exits a bounded map 
(j)- 1 : H —> C (Hl 3 _ 5 ; H) such that for any x £ H, 

4>(x) ■ </> _1 (a:) • h = h for any h £ H 3 _ 5 . 

Moreover f £ V. 

The method to extract (Nk)k is based on the investigation of the properties of the 
Kolmogorov equation associated to o> perturbed by a very irregular potential. 

It follows that (Px 0 )a;oeH 2 is a weak Markov family, which means that for any 
x 0 £ H 2 

(1.11) IP X Q (X*(t) £ H 2 ) = 1 for any t > 0. 
and that, for any t\ < ■ ■ ■ < t n , t > 0 and any if £ B b ( H 2 ; R) 

(1.12) E* 0 (^pC.(i + t n ))| X*(H),... ,X*(t n )) = Vt^X^tn)), 
where 

(Vtif) (x 0 )=E Xo (ip(X*m- 

Note that OH) was known only for a stationary solution (see 1X3). 

Remark 1.11. Assume that Hypothesis o holds. If we strengthen ns into 

K2 = ^l <0 °’ 

for some 5 > 0, then Hypothesis Oil holds. 

212 





Exponential mixing for the 3 D stochastic Navier-Stokes equations 


Hence, we immediately deduce the following corollary from Theorem ll. 8 l 

Corollary 1.12. Assume that Hypothesis E3 and, \T~M hold. Then there exits a 
unique invariant measure /.i for (7 7 t) tgK + and C, 7 > 0 such that for any A £ 'P(H 2 ) 


(1.13) ||P t *A - p\\ var < Ce-* ( 1 + f M 2 X(dx)) , 

\ Jm 2 J 

provided ||/|| 2 < <5° and where is the total variation norm associated to the 

space H 2 . 

Remark 1.13 (Uniqueness of the invariant measure p). Assume that Hypothesis 
rrm holds. Let P Xo and P(, o be two weak solutions of o which are limit in 
distribution of solutions of O- Then we build (Vt)t an d (Pt)t as above associated 
to IP^q and vespectxvely. It follows that theve exzsts /x and [x such that (trial) 
and ( 11 . 101 ) hold for P Xo ,/z) and ({V' t )t, P(. Q , p'). Although we have uniqueness 

of the invariant measures p and p! associated to (Vt)t and (' P[)t> we do not know 
if p and p! are equal. 


1.4. Coupling methods. 

The proof of Theorem Ol is based on coupling arguments. We now recall some 
basic results about coupling. Moreover, in order to explain the coupling method in 
the case of non degenerate noise, we briefly give the proof of exponential mixing 
for equation O- 

Let (Ai, A 2 ) be two distributions on a polish space (E,£) and let (12, T, P) be a 
probability space and let (Zi,Z 2 ) be two random variables (12, T) —» (E,£). We 
say that (Zi,Z 2 ) is a coupling of (Ai,A 2 ) if A i = V(Zi) for i = 1,2. We have 
denoted by T>{Zi) the law of the random variable Z,. 

Next result is fundamental in the coupling methods, the proof is given for in¬ 
stance in the Appendix of ESI. 

Lemma 1.14. Let (Ai, A 2 ) be two probability measures on (E,£). Then 

IIAi - A 2 ||,, Qr = minP(Zi ^ Z 2 ). 

The minimum is taken over all couplings (Zi,Z 2 ) of (Ai, A 2 ). There exists a cou¬ 
pling which reaches the minimum value. It is called a maximal coupling. 


Let us first consider the case of the solutions of O- Assume that Hypothesis 
o holds. Let N £ 1A and (xq,Xq) £ R 2 . Combining arguments from ESI, EH, it 
can be shown that there exists a decreasing function pn(-) > 0 such that 


(1.14) 


0 P l)* 5 xl - ('PlY 6 xl < 1 -PN (Kl + l^ol) 

u u 11 rtnr- 


Applying Lennna ri.141 we build a maximal coupling (Z\, Z 2 ) = ( Z\(xq , Xq), Z 2 (xq, Xq 
of ((Pf)* 5 X 1 , (V^Y S x 2 ). It follows 

(1.15) P(Z 1 = Z 2 )>p N (\x 1 o \ + \xl\)>0. 


Let (IT, W) be a a couple of independent cylindrical Wiener processes and <5 > 0. 
We denote by Xn{-,xq) and Xjy(-,xo) the solutions of (11.311 associated to W 
and W, respectively. Now we build a couple of random variables (Vi,V 2 ) = 
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(Vi(xq,Xq), V 2 (xq,Xq)) on PnH as follows 
(X N (-, x 0 ), X N (•, x 0 )) 

(Zi(xq, Xq), Z 2 (xq, Xq)) 

(X N (-,xl),X N (-,x%)) 


(1.16) (V 1 ,V 2 )={ 


if xl = xl = x 0 , 

if (xl,xl) e B ff (0,(5)\{a;J = x%}, 
else, 


where Bh(0,5) is the ball of H x H with radius 6. 

Then x%), V 2 (xl, z§)) is a coupling of ((Vi)*5 x i, (P^)* S x 2 ). It can be 

shown that it depends measurably on (xI,Xq). We then build a coupling (X 1 , X 2 ) 
of (V(Xn(-,Xq)),V(Xn(-,Xq))) by induction on N. We first set A'*(0) = x l 0 for 
i = 1,2. Then, assuming that we have built (X^X 2 ) on {0,1,..., k}, we take 
(Vl,V 2 ) as above independent of (X^X 2 ) and set 

X\k + 1) = V)(X 1 (fc),X 2 (fc)) for i = 1,2. 


Taking into account <0>, it is easily shown that the time of return of (X^X 2 ) 
in B(0,4:(c/ni)Bi) admits an exponential moment. We choose i5 = A(c/n\)B\. It 
follows from 11.151) . II. Kill that, (X' 1 (n), X' 2 (n)) G B(0, S) implies that the proba¬ 
bility of having (X 1 , X 2 ) coupled (i.e. equal) at time n + 1 is bounded below by 
Pat(2i 5) > 0. Finally, remark that if (X 1 , A' 2 ) are coupled at time n + 1, then they 
remain coupled for any time after. Combining these three properties and using the 
fact that (X 1 (n), X 2 (n)) nS N is a discrete strong Markov process, it is easily shown 
that 

(1.17) P (X 1 (n) ^ X 2 (n)) < C N e~™ n (l + + |;rg| 2 ) , 

with 7 jv > 0 . 

Recall that (X' 1 ,X 2 ) is a coupling of (P(XW(-,^o)),T>(X at(-, a; 2 ))) on N. It 
follows that (X' 1 (n),X' 2 (n)) is a coupling of ({Vn)*8 x i, ('P ? ^ r )*<5 x 2). Combining 
Lemma im and 03. we obtain, for n E N, 


{K)*8 xl - (VZ)*6 X 1 < C N e~™ n (l + \xl\ 2 + \x 2 0 \ 2 ) 

u u var \ / 


Setting n = [_ij and integrating ( Xq,x ^) over ((V^_ n )*X) g) Hn where is an 
invariant measure, it follows that, for any A £ P(PnH) 1 


(1.18) 


C Pt) ^ ~ 


< C N e~’ 1Nt 

var 



This result is useless when considering equation o since the constants Cn, 'Jn 
strongly depend on N. If one tries to apply directly the above arguments to the 
infinite dimensional equation o, one faces several difficulties. First it is not 
known whether P Xo is Markov. We only know that, as explained in section 1.3, a 
Markov transition semi-group can be constructed. This is a major difficulty since 
this property is implicitely used in numerous places above. Another strong problem 
is that Girsanov transform is used in order to obtain irm Contrary to the two 
dimensional case, no Foias-Prodi estimate is available for the three dimensional 
Navier-Stokes equations and the Girsanov transform should be done in the infinite 
dimensional equation. This seems impossible. We will show that we are able to 
prove an analogous result to itm by a completely different argument. However, 
this will hold only for small initial data in H 2 . Another problem will occur since it 
is not known whether solutions starting in H 2 remain in H 2 . 
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We remedy the lack of Markov property by working only on Galerkin approx¬ 
imations and prove that (P 8 l) holds with constants uniform in N. As already 
mentioned, we prove that is true for Xq, Xq in a small ball of H 2 and uni¬ 

formly in N. Then, following the above argument, it remains to prove that the time 
of return in this small ball admits an exponential moment. Note that the smallness 
assumption on / is used at this step. In the following sections, we prove 


Proposition 1.15. Assume that Hypothesis o holds. Then there exist 6° = 
5 °{Bq, D, e, v), C = C((j>, D , £,i/)> 0 and 7 = 7 (</>, D, e, v) > 0 such that if \\f\\1 < 
6° holds, then, for any N G N, there exists a unique invariant measure hn for 
('Pf , ) tm+ - Moreover, for any A G V(PnH) 


(1.19) 


(V t N )*\-» N 


< Ce~^ 



We now explain why this result implies Theorem 1 1.81 

Let A G P(H) and X\ be a random variable on H whose law is A and which is 
independant of W. Since ||• ||^, ar . is the dual norm of l-^, then (TTTfll) implies that 


( 1 . 20 ) 


K(g{X N {t,X x )))- 


' PnH 


g{x) n N (dx) 


<C\g \ c 


1 + / \x\ A (dx) 
Jh 


for any g G UCb(M s ) for s < 0. 

Assume that, for a subsequence (N' k )k, Xjsi(t,X\) converges in distribution in 
H s to the law X*(f) under the weak solution of 11 .111 . Recall that the family 
00 jv is tight. Hence, for a subsequence (Nk)k of {N k )k, P^ Wfc converges to P^ a 
weak stationary solution of o> with initial law /i. Taking the limit, Ol follows 
from 11.2011 . which yields Theorem II.81 


2. Coupling of solutions starting from small initial data 


The aim of this section is to establish the following result. A result analogous to 
(tTTKt but uniform in N. 


Proposition 2 . 1 . Assume that Hypothesis E3 holds and that f G H. Then 
there exist ( T,S ) G (0, l ) 2 such that, for any N G N, there exists a coupling 
(Zi(xq,Xq), Z 2 ^Xq,Xq)) of{{V^)*S x i, [V^)*5 X 2) which measurably depends on (xq,Xq 
H 2 and which verifies 

( 2 . 1 ) F(Z 1 (xl,xl) = Z 2 (xl,xl)) > ^ 

provided 

(2.2) IKH 2 v IKHa < S. 


Assume that Hypothesis 11.41 holds and that / G H. Let T G (0,1). Ap¬ 
plying Lemma irui we build (Zi(xq, Xq), ^ 2 (^ 0 , £§)) as the maximal coupling of 
(Pj’Sj.i , PjiS x 2 ). Measurable dependance follows from a slight extension of Lemma 
1.17 (see remark A.T). 

In order to establish Proposition l 2 H it is sufficient to prove that there exists 
c(Bq,D) not depending on T G (0,1) and on TV G N such that 


(v»y 5xl -(v?y5 xl 
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< c{B 0 ,D)Vt, 


(2.3) 


var 
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provided 

(2-4) IKHaVHH 1<B 0 T 3 . 

Then it suffices to choose T < l/(4c(B 0 , D)) 2 and S = B 0 T 3 . 

Since IHI^^ is the dual norm of |- EH is equivalent to 

(2.5) |E (g(X N (T,xl)) — g(X N (T,xl))) | < 8 c(B 0 , D)Vt. 

for any g £ UCb(PwH). 

It follows from the density of C^(PnH) C UCb(PNH) that, in order to establish 
Proposition EH it is sufficient to prove that 12.511 holds for any TV £ N, T £ (0,1) 
and g £ CI(PnH) provided 12.411 holds. 

The proof of 12.511 under this condition is splitted into the next three subsections. 

2.1. A priori estimate. 

For any process X, we define the Hi-energy of X at time t by 

Ef(t) = \\X(t)\\ 2 + f \\X{s)\\lds. 

Jo 

Now we establish the following result which will be useful in the proof of 12.51 

Lemma 2.2. Assume that Hypothesis E3 holds and that f £ H. There exist 
Ko = Kq(D) and c = c(D) such that for any T < 1 and any N £ N, we have 

provided ||xo|| 2 < BfT. 

Let Xn = Xn(-,xq). Ito Formula on ||Xjv|| 2 gives 

(2.6) d ||ATjv|| 2 + 2 ||X/vlj 2 dt = dM Hl + I^dt + ||-Pjv</(Xjv)||£ 2 ( i?;Hl ) dt + Ifdt , 

where 

j J Hl = -2 {AX n , B{X n )) , I f = 2 (AX n , f ), 

| M Ul (t) = 2 f* ( AX N (s ), <t>(X N (s))dW(s)). 

Combining a Holder inequality, a Agmon inequality and a arithmetico-geometric 
inequality gives 

(2.7) / Hl <2 ||X*|| 2 |X*U|X*|| <c\\X N \\l p^ll* <±\\X N \\ 2 2 + c\\X N \\ 6 . 
Similarly, using Poincare inequality and Hypothesis 11.41 

(2.8) I f <\ ||Xjv II 2 + c |/| 2 < \ ||X^|| 2 + cB 0 . 

We deduce from EU- E3b (ESI. Hypothesis O and Poincare inequality that 

(2.9) d ||X/v|| 2 + IIXjv-H \dt < dM ai + cB 0 dt + c HXjvf (llX^f - 4X 0 2 ) dt, 
where 



( 2 . 10 ) 
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Setting 


a Hl = inf [t £ (0, T) | ||Wv(t)|| 2 > 2A' 0 j, 


we infer from ||a:o|| 2 < BqT that for any t £ (0,cthi) 

(2.11) Ef N (t)<cB 0 T + M Ml (t). 

We deduce from Hvpothesis ll ,4l and from Poincare inequality that A is bounded 
in £(Hi;IHi) by cB 0 . It follows that for any t £ (0,aHi) 

(M Hl )(i)=4 / \\P N <j>(X N (s))*AX N (s)fdt <cB 0 [ \\X N (s)fds<2cK 0 B 0 T. 
Jo Jo 

Hence a Burkholder-Davis-Gundy inequality gives 

E ( sup M Hl ) < cE^(M Hl ) (cthJ < c\JK 0 B 0 T < c(K 0 + B 0 )VT. 

\(0.ct H i ) J 

It follows from (12.111) and T < 1 that 

E ( sup Ef N ) < c(B 0 + K 0 )Vf, 

\( 0 , ct Hi ) ) 

which yields, by a Chebyshev inequality, 


( sup <c(l + ^)vf. 

\(0,a Wl ) ) V K 0 J 


Now, since sup( 0crHi ) A®^ < Kq implies cth! = T , we deduce Lemma 12.21 


2.2. Estimate of the derivative of Ij\r. 

\2 


Let TV £ N and ( Xq , h) £ (Hb) . We are concerned with the following equation 

dijN + Arj N dt + P n B(X n ,t]n) dt = Pn^'(Xn) • Vn) d,W, 
r) N (s, s,x 0 ) ■ h = P N h, 


( 2 . 12 ) 


where B(X N ,rj N ) = B(X N ,r ) N ) + B(r] N ,X N ), X N = X N (-,x 0 ) and i] N (t) = 
i]N{t, s , Xq) ■ h for t > s. 

Existence and uniqueness of the solutions of itm are easily shown. Moreover 
if g £ CI(PnH), then, for any t > 0, we have 

(2.13) (V (V t N g) (x 0 ),h) = E (Xg(X N (t, x 0 )), m(t, 0, *o) • h ). 

For any process X , we set 

rt 


(2.14) 


a(X) = inf it £ (0, T) 


f ||X(s)|| 2 ds> K 0 + l 
Jo 


where Kq is defined in Lemma o We establish the following result. 

Lemma 2.3. Assume that Hypothesis \l . f\ holds and that f £ H. Then there exists 
c = c(I3o, D) such that for any N £ N, T < 1 and (xo, h) £ (H 2) 2 

na{X jv (-, x 0 )) 

E / 


||77vr(t,0, x 0 ) ■ h\\l dt < c" L " 2 


2 • 
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For a better readability, we set ijN(t) = VN(t,0, xq) ■ h and a = <j(Xn(-,xo)). Ito 
Formula on ||r7Ar(t )||2 gives 

(2.15) d\\r) N \\l + 2 \\r] N \\ldt = dM VN + I m dt+\\P N (<j)'(X N ) ■ Vn)\\c 2{u . B2) dt , 

where 

I Mwd) = 2 /o { A2 VN,{PN<t>'{X N )-ri N )dW) ds, 

\ I m = -2 (a*t, n ,A$B( X N , m )y 

It follows from Holder inequalities, Sobolev Embedding and a arithmetico-geometric 
inequality 

< c 1177/v11 3 ||t7jv || 2 ll-^Jvl^ ^ + c II 7 ?Jv ||2 II^Jvllz • 

Hence, we deduce from I2.15K and Hypothesis 11.41 


d \\Vn\\1 + \\Vn\\1 dt < dM m + c \\rj N \\l \\X N \\l + B 0 \\mWl dt. 

Integrating and taking the expectation, we obtain 

(2.16) E^S{a,0)\\r] N (a)\\ 2 2 + £((T,t)\\r] N {t)\\l dtj < \\h\\l, 

where 

£(t,s) = H Xjv i r )ll 2 dr . 

Applying the definition of ct, we deduce 

(2.17) E f \\ri N (t)\\ldt<\\h\\lexp(c(K 0 + l) + B 0 T), 

J 0 

which yields Lemma EH 


2.3. Proof of (12.511 . 

Let ip £ C°° (R; [0,1]) such that 

ip = 0 on (Lfo + 1) 00 ) an d ip = 1 on (—oo,A' 0 ). 
For any process A', we set 

Remark that 

(2.18) |E (g(X N (T, ®g)) - g(X N (T, *£))) \<I 0 + ^ (h + I 2 ), 
where 


do — 
h = 


E ( g{X N (T,xl))tpx N (;xl ) - g{X N (T,x\))ipx N (-,xl)) 
1 (/ 0 T || Aat(s, a?o)||2 ds > K o) ■ 


For any 6 £ [1, 2], we set 


%o = (2 - 0)x J + (6- l)xl, X e = X N (-, xl), 


m{t) = VN(t, 0,xg), 


a g = cr(X g ). 
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Recall that a was defined in 1 EH) . For a better readability, the dependance on N 
has been omitted. Setting 

h Xq Xq , 

we have 


(2.19) h ~j l |Je|d0 Je = (VE(g(X e (T))i> Xe ),h). 

To bound Jg , we apply a truncated Bismut-Elworthy-Li formula (See appendix fXl) 


( 2 . 20 ) 

where 

p 

J e,i 

< J'e, 2 


E {g(Xg(T))^ Xg ■ Ve (t) ■ h,dW(t ))) , 

E (g{X 0 (T)W Xe /; e (1 - 1) (dX fl (t), A( Ve (t) ■ h)) dt) , 
*!>' (f 0 T \\Xe(s)\\lds) . 


It follows from Holder inequality that 


\ J h\ ^ l^loo l^'loo -y 

and from Hypothesis 11.41 that 


E 


\\Xg(t)\\ 2 dtJ E / II I)g{t) • h\\ 2 dt. 


I T' I 
\ j ba\ 


<lfflc 



II m(t) ■ Hi dt. 


Hence for any T < 1 


(2.21) \J g \ < c(B 0l D) 1.9^ ^E jH II ne{t) • dt. 

Combining and Lemma roi we obtain 

W<c(B 0 ,D)|s| oo ll^ > 

which yields, by m and OTT5I . 

I 0 <c(B 0 ,D) | g^VT. 

Since BqT 3 < BqT, we can apply Lemma 12 . 21 to control Ji + I 2 in iTHsli if (IQ) 
holds. Hence E3 follows provided (IOi holds, which yields Proposition 12. II 


3. Time of return in a small ball of H 2 

Assume that Hypothesis 11.41 holds. Let N £ N and T,S,Zi,Z 2 be as in Propo¬ 
sition l2Tl Let (IT, VP) be a couple of independant cylindrical Wiener processes 
on H. We denote by Xn(-,xo) and Abv(-,£o) the solutions of HI.311 associated 
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to W and W, respectively. We build a couple of random variables (14,14) = 
(Vi(xq, x%), V 2 (xl, Xq)) on PnH as follows 


(3.1) (V 1 ,V 2 )={ 


(Xjv(-, xo ), X N (-, x 0 )) 
(Z 1 (x^,x§),Z 2 (x^,x§)) 
(X N (-, Xq), Xn(-, Xq)) 


if xl = Xq = x 0 , 

if (xj, xl) e Bm 2 (0, 6)\{xl = x§}, 
else, 


We then build (X^X 2 ) by induction on TN. Indeed, we first set X l (0) = x z 0 for 
i = 1,2. Then, assuming that we have built (X^X 2 ) on {0, T, 2T, ..., nT}, we 
take (Vi, 14) as above independent of (X^X 2 ) and we set 


X*((n + 1)T) = Vi(X 1 (nT), X 2 (nT)) for i = 1,2. 

It follows that (X'-^X 2 ) is a discrete strong Markov process and a coupling of 
(2?(Xjv(-, Xq)), 2?(Xjv(-, Xq))) on TN. Moreover, if (X^X 2 ) are coupled at time 
nT, then they remain coupled for any time after. 

We set 


(3.2) 


inf jt £ TN\{0} ||X 1 (t)||^v||X 2 (t)| 



The aim of this section is to establish the following result. 


Proposition 3.1. Assume that Hypothesis [T^\ holds. There exist 6 3 = 6 3 (Bq, D,e,5), 
a = a{4>,D,£,8) > 0 and K" = K” (</>, D, e, <5) such that for any (xJ,Xq) £ H x H 
and any N £ N 

E (e aT ) < K" (l+|x(,| 2 + |x2| 2 ), 

provided \\f\\1 < 5 3 . 


The result is based on the fact that, in the absence of noise and forcing term, 
all solutions go to zero exponentially fast in H. A similar idea is used for the two- 
dimensional Navier-Stokes equations in ES|. The proof is based on the following 
four Lemmas. The first one allows to control the probability that the contribution of 
the noise is small. Its proof strongly uses the assumption that the noise is diagonal 
in the eigenbasis of A. As already mentioned, in the additive case, the proof is easy 
and does not need this assumption. 


Lemma 3.2. Assume that Hypothesis O holds. For any t,M > 0, there exists 
Po(t, M) = po(t, M, e, (|(/'n| 00 )n, D) > 0 such that for any adapted process X 


where 


¥ I sup \\Z\\ 2 2 <m) 
\(° ■*) / 


>p 0 (t,M), 



e - A{t - s U(X(s))dW(s). 


It is proved in section 3.1. 

Then, using this estimate and the smallness assumption on the forcing term, we 
estimate the moment of the first return time in a small ball in H. 

Let <5 3 > 0. We set 


tl 2 = r A inf |t € TN* 


X^fvIX 2 ^)! > <5 3 j . 
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Lemma 3.3. Assume that Hypothesis \1.4\ holds. Then, for any 83 > 0, there exist 
03 ( 83 ), 03 ( 83 ) and 73 ( 83 ) such that for any (xq,Xq) £ (H 2) 2 

E(e^)<C 3 (l + \x 1 0 \ 2 + \x 2 0 \ 2 ), 

provided 

I/I < C' z . 

The proof is postponed to section ESI 

Then, we need to get a finer estimate in order to control the time necessary 
to enter a ball in stronger topologies. To prove the two next lemmas, we use an 
argument similar to one used in the determinist theory (see BS], chapter 7). 


Lemma 3.4. Assume that Hypothesis E3 holds. Then, for any 8 4 , there exist 
Pa(8a) > 0, 64 ( 154 ) > 0 and Ri( 8 i) > 0 such that for any Xq verifying |cco | 2 < R4, 
we have for any T < 1 

p(||X J v(7>o)|| 2 <<54) >p 4 , 


provided 


I/I < ci. 


The proof is postponed to section Ena 


Lemma 3.5. Assume that Hypothesis E3 holds. Then, for any S 3 , there exist 
P 5 ^ 5 ) > 0, 03 ( 85 ) > 0 and i?s(^ 5 ) > 0 such that for any Xq verifying ||aro|| 2 < R 5 
and for any T < 1 

P (\\X n (T, * 0 ) 1)2 < ^ 5 ) > P5- 

provided 

\\f\\ £ < C h . 

The proof is postponed to section EH1 


Proof of Provosition ld.ll We set 

85 = 8 , 84 = R5(85), S 3 = R4(84:), P4=P4(80), P5 = P5(^), Pi = (P4Pdf , 

and 

8 3 = 0 ( 3 ( 83 ) ACM A C' 5 ( 8 5 ). 

By the definition of t l 2 , we have 

l ^ 1 ^,)! 2 v | js : 2 (^ 2 )| 2 < ^ 4 (^ 4 ). 

We distinguish three cases. 

The first case is ||A' 1 (r i 2)|| 2 V ||A 2 (r i 2)|| 2 < 8 , which obviously yields 

(3.3) P ^ min ^ max || A' 1 (t L 2 + fcT )|| 2 < (5 | (X 2 (t L 2 ), X 2 (t L 2 )) ^ >p\. 

We now treat the case *0 = A 1 (r i 2 ) = A 2 (tl 2 ) with \\xq\\\ > <5- Combining 
Lemma l3.4l and Lemma, 13.51 we deduce from the weak Markov property of Xn that 

P (11^(27,^)112 <*) >P5P4, 
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provided |.To | 2 < Ra- Recall that, in that case, X 1 (t L 2 +2T) = X 2 (t L 2 +2T). Hence, 
since the law of X l {r L 2 + 2 T) conditioned by {X\t L 2),X 2 (tl 2 )) is D(Xn( 2 T, xo)), 
it follows 

P (™i2 11** fa 2 + 2T )\\2 - S I ( xl ( T i 2 )> A ' 2 ( T i 2 )) ^ > PAP5 > Pi, 
and then (TOll 

The last case is X l {r L 2) ^ X 2 (tl 2 ) and ||( 7 -^ 2 )||^ V ||X 2 (r L 2)|| 2 > S. In that 
case, (X 1 (tl2 +T),X 2 (t L 2 +T)) conditioned by (X 1 (t L 2),X 2 {tl 2 )) are indepen¬ 
dent. Hence, since the law of X l (r L 2 + T ) conditioned by (X 1 (r i 2 ), X 2 {tl 2 )) = 
(xq,Xq) is V(Xn(T,x 1 0 )), it follows from Lemma mi that 

P (max ||X l (r i 2 +T )|| 2 < | (X 1 (r L2 ), X 2 (tl 2 )) ^ > p\. 

Then, we distinguish the three cases (||X*(t £, 2 + T)|| )j = i i2 in the small ball of H 2 , 
equal or different and we deduce from Lemma Id. 51 bv the same method 

P f min max ||X J (r i2 + kT) || 2 < S | ( 2 Y 1 (r i 2 + T),X 2 (t L2 + T)) j > p\, 

\ k — 1,2 i — 1,2 / 

provided 

max II X % (tl,2 + T) II << 54 . 

2=1 2 1 

Combining the two previous inequalities, we deduce mi for the latter case. We 
have thus proved that mi is true almost surely. 

Integrating mi . we obtain 

(3.4) P ( min max 11X 1 (77,2 + fcT )|| 2 < S ) > m. 

\k=0,...,2i=l,2 11 112 J 

Combining Lemma 13.31 and m, we conclude. 

3.1. Probability of having a small noise. 

We now establish Lemma mi 

We deduce from Holder inequality and from J2 n P'n 2 < 00 that Hypothesis 11.41 
implies the following fact: for any e 0 £ ( 0 , e), there exists a € ( 0 , 1 ), a family (</>„)„ 
of measurable maps H —> K. and a family (&,), of positive numbers such that 

J 

\ sup xGH \<i>n{x)\ <1, B* = J 2 n /^ +e ° ( b n f {1 ~ a) < 00. 

For simplicity we restrict our attention to the case t = 1. The generalization is 
easy. 

Remark that 

Z(t) = y ^b n Z n (t)e n , 

n 

where 

Z n (t) = [ e~ llnl ' t ~ s ' > 4> n (X(s))dW n (s), where W = Y]W n e n . 

Jo 
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It follows from \\Z\\ 2 = Yhn^nTn \ M^Z n \ 2 and from 13.511 that 

(3.6) P ( sup \\Z\\l < B*M ] > P ( sup \y^Z n \ 2 < MpQ ( b n )~ 2a , V n ] . 


,( 0 , 1 ) 


,( 0 ,!) 


Setting 


W„{t) = yfJI^Wn , 

we obtain (Wf) n a family of independent brownian motions. Moreover we have 

\fP : nZ T i (t) = Z n (fl n t), 

where 

Z' n (t) = j o e-«- S ^ n (s)dW^( S ), ^ n (s) = K (* 

Hence, it follows from cut that 


(3.7) 


sup ||Z||* < B*M > P sup \Z' n \ 2 < MpQ ( b n )~ 2a , V1 


,(0,1) 


l(0,/in) 


Let Wl i = W' n (i + ■) - W^i) on (0,1). We set 


{ 0 if t < 0, 

f 0 1M e s V’n (i + s)dW ’ n i (s) if t > 0. 

Remark that 

OO 

i =1 

which yields for any q G N 

(3.8) sup \\Z’ n \\ 2 < (— S — ] max sup | M n ,i\. 

(o, 9 ) \e — iy t=o,...,9-i ( 0 , 1 ) 

Remark that [W' nk ) n ,k is a family of independant brownian motions on (0,1). 
It follows that ( M nt k)n,k are martingales verifying (M n ^, M n < ik ') = 0 if (n, k) ^ 
{n',k'). 

Hence, combining a Theorem by Darnbis, Dubins and Schwartz (Theorem 4.6 page 
174 of I2IJ) and a Theorem by Knight (Theorem 4.13 page 179 of J233J), we obtain 
a family {B n ^) n ,k of independent brownian motions verifying 

(3.9) M„ ik (t) = Bn,k((M n> k) (i)). 

Remark 3.6. In the two previous Theorem (Theorem 4-6 page 174 and Theo¬ 
rem 4-13 page 179 of m). it is assumed that P (< M > ( 00 ) < 00 ) = 0. How¬ 
ever, as explained in Problem 4 -7 of m, the proof is easily adapted to the case 
P (< M > ( 00 ) < 00 ) > 0. 

Remarking that for any t £ (0,1) 
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we deduce from (|3.8|) and 13.911 that for any q £ N* 

sup \\Z ' n \\ 2 < ( ——- ) . max sup \B n ^\. 

(0 ,q) V e ~ 1/ *=0,...,g— 1 (o,l) 

Hence it follows from ED that 


sup \\Z \\ 2 <M I > P I sup \B n ^\ z < cM /4° (b n ) 2 “, Vn, V*</i„ + 1 I , 


do. 1 ) 


,(o,i) 


where c=(^i) 

We deduce from the independence of {B n ^) n ,k that 


(3-10) P (jup ||Z||“ < M) > [] (p ( cM ( b n )~ 2a ) 
where 


raSN* 


Mn + 1 


P(d 0 ) = P sup l-Bqif < d 0 . 

V(o.D ) 

Recall there exists a family (c p ) p such that 

E ( sup |Pi,i| 2p ] < c p . 

V(o,D J 

It follows from Chebyshev inequality and from 1 — x > e~ ex for any x < e _1 that 

1 

for any do < d p = (e _1 c p ) p 

P(d 0 ) > 1 - c p do P > e~ ec * d * V . 

Applying 13.101) . we obtain for any p > 0 


(3.11) P ( sup \\Z\\\ < M 1 > C p (M)ex p I — 


,(o,i) 


Mp 


E 


n>N(p,M) 


l^n + 1 


£o P I n 


bl ap 


where 


N{p, M) = sup | n e N\{0} M n £ ° (b n ) 2 “ < d p |, 

c p {M) = n„< N(PiM) [p (cMvtf ( b n y 2a Yj . 

Choosing p sufficiently high, we deduce from HO that 

\ f'ri J 
n x 7 

which yields, by (THT1) . that for any M > 0 and for p sufficiently high 


(3.12) 


sup ||^||2 < M J > C P (M) exp ( — 

,(o,i) / 


Mp 


Remark that for any p , £q we have N(p, M) < oo. Moreover, it is well-known that 
for any do > 0, P(do) > 0, which yields C P (M) > 0 and then Lemma 1301 
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3.2. Proof of Lemma EH 

For simplicity in the redaction, we restrict our attention to the case f = 0. The 
generalisation is easy. 

Recall (11.511 

E |X/v(f )| 2 < e _Ml * |cco | 2 + —Bq. 

Mi 

Since (X^X 2 ) is a coupling of (V(Xn(-,Xq)),V(Xn(-,Xq))) on TN, we obtain 

(3.13) E (|A’ 1 (nT)| 2 + |X 2 (nT)| 2 ) < e~^ nT (jxjj 2 + |zq| 2 ) + 2^-R 0 - 

Since (X 1 , X' 2 ) is a strong Markov process, it can be deduced that there exist Cq 
and 76 such that for any Xq G H 

(3.14) E (e^ 2 ) < C 6 (l + |*J | 2 + |ag| 2 ) , 


where 


t ' l2 = inf {f g TN\{0} |^T 1 (t)| 2 + \X 2 {t )\ 2 > 4cR 0 } 


Taking into account a standard argument gives that, in order to establish 

Lemma 1301 it is sufficient to prove that there exist (py,Ty) such that 


(3.15) F a; 0 )| 2 < J 3 ) >Pr(<M)> 0, 

provided N G N, t > ^( 63 ) and |cco | 2 < 4c.Bo- 
We set 

Z{t)= [ e- A ^- s U(X N {s))dW{s), Y n =X n -P n Z, M = sup||Z|| 2 . 
Jo (0 ,t) 


Assume that there exist My( 83 ) > 0 and Ty(5s) such that 

(3.16) M < My(S 3 ) implies |Yjv(i)| 2 < 

provided t > Ty(Ss) and |xo|“ < 4c.Bo- Then 13.1511 results from Lemma 1301 with 

M = min |m 7 (<5 3 ), ^ j . 

We now prove EE3. Remark that 

(3.17) jY n + AY n + P n B(Y n + P N Z) = 0. 

Taking the scalar product of (I3H71) with Y n , it follows that 

(3.18) j t IFjvI 2 + 2 ||y/v|| 2 = -2 (Y n , B(Y n + P N Z)). 

Recalling that (B(y,x),x) = 0, we obtain 


-2(Fjv, B{Y n + P N Z)) = —2(Y n , (Y n , X)P n Z ) - 2 (Y N , B{P N Z)). 


We deduce from Holder inequalities and Sobolev embedding that 

~(z,(x,V)y) <c\\z\\ ||x|| ||y|| . 

Hence it follows from (13081) that 

j t |Yat| 2 + 2 ||F W || 2 < c||Z|| 2 \\Y n \\ + c\\Z\\ HYvII 2 , 
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which yields, by an arithmetico-geometric inequality, 

j t IFjvI 2 + 2 ||1jv|| 2 < cMi ||y}v|| 2 + cAf*. 

It follows that M < 4j- implies 

(3.19) ^\Y N \ 2 + \\Y N \\ 2 <cMi on (0 ,t). 

at 

Integrating, we deduce from |xo| 2 < 4ci?o that 

\Y N (t)\ 2 < 4ce _/il *Bo +c ' 

Choosing t sufficiently large and M sufficiently small we obtain (13. Kill which yields 
(13.1511 and then Lemma roi 

Remark 3.7. In order to avoid a lengthy proof , we have not splitted the arguments 
in several cases as in the proof of Proposition \d.l\ The reader can complete the 
details. 


3.3. Proof of Lemma rr~n 

We use the decomposition Xn = Yn + PnZ defined in section IQ and set 

M = sup \\Z\\l. 

(0,T) 

Integrating (EHTl . we obtain for M satisfying the same assumption M < 

Jo ^ Yn ^ 2 dt ~ \ ^°! 2 + cM ^’ 

which yields, by a Chebyshev inequality, 


(3.20) 


A ie(0,T) 


\Y N (t)\\ 2 <^\x 0 \ 2 + 2cM-* J >|, 


where A denotes the Lebesgue measure on (0,T). 
Setting 


r H i = inf t e (0, T) 


\Y N (t)\\ 2 <^\x 0 \ 2 + 2cMi 


we deduce from eoto and the continuity of Yn that 
(3.21) \\Y N ( mi )f<^\x 0 \ 2 + 2cMi. 

Taking the scalar product of 2AY and EH7I) . we obtain 


(3.22) j t || Vac || 2 + 2 \\Y N f 2 = -2(AY N , B{Y n + P N Z )). 

It follows from Holder inequalities, Sobolev Embeddings and Agrnon inequality that 

-2 (Ay,B(x,z)) < c\\y \\ 2 ||z|| 5 ||z||f ||x|| , 
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where B{x,y) = (x,'V)y + ( y,V)x. Hence, we obtain by applying arithmetico- 
geometric inequalities 


-2 (AY n ,B(Y n )) 

< 

c\\Yn\\I \\Y N p 

< 

l\\Y N \\ 2 2 + c\\Y N \\ 6 

-2 (AY n ,B(P n Z)) 

< 

c \\Y N \\ 2 \\Z\\i \\Z\\I 

< 

\\\Y N \\ 2 2 + c\\Z\\l, 

—2(AY n , B(Y n , P n Z)) 

< 

c\\Y N \\l \\Y N p \\Z\\ 

< 

c\\Z\\ \\y n \\ 2 2 


Remarking that B(Yn + PnZ) = B(Yn) + B(Yn, PnZ) + B(PnZ ), it follows from 
(ET221) that M < implies 

(3.23) | Hl^vll 2 + \\Y N f 2 < c\\Y n \\ 2 (||l^v|| 4 - 4 K 2 ) + cM 2 , 

where Kq is defined in iTHUl) . Let us set 

a Hl = inf [t e (r Hl ,T) | ||^v(t)|| 2 > 2K 0 } , 

and remark that on (th! , <tm 1 ), we have 

d 


(3.24) 


dt 


IIL/vIT + \\Y n \\1 < cM 2 . 


Integrating, we obtain that 

(3.25) ||Liv( ( 7H 1 )|| 2 + / \\Y N (t)\\ldt<\\Y N (T Bl )\\ 2 +cM 2 . 

Combining iEhti and (ITT251) . we obtain that, for M and |a;o| 2 sufficiently small, 

II^hJII 2 < 

which yields <thi = T. It follows that 

(3.26) IIAbvCOH 2 < <S 4 , 

provided M and |xo| sufficiently small. It remains to use Lemma l3.2l to get Lemma 

IQ 

3.4. Proof of Lemma EH 

It follows from <EZ1 that 

[ Il^v(0ll2 dt ^ ll^oll 2 + cM 2 , 

Jo 

provided M < A and ||xo|| 2 + cM 2 < K 0 - 

Applying the same argument as in the previous subsection, it is easy to deduce that 
there exists a stopping times th 2 G (0, T) such that 

(3.27) ' '" 2 ' 2 


\Y N (m 2 )\\ 2 2 <- (\\x 0 f+cM 2 ) 


provided M and ||xo|| are sufficiently small. 

Taking the scalar product of 13.1 7ll and 2A 2 Yjv, we obtain 

d ||Ljv|| 2 + 2 HYaHI 2 = —2 (A*Y N ,AlB(y N + P N Z )) . 


(3.28) 


dt 
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Applying Holder inequality, Sobolev Embeddings H 2 C L°° and Hi C L 4 and 
arithmetico-geometric inequality, we obtain 

-2 (A?y,A?B{x,y)} < c\\y\\ 3 ||x|| 2 ||y|| 2 < i \\yf 3 + c (||:r|| 2 + \\y\\fj . 

Hence we deduce from (13.281) and from B(Yn + PnZ) = B{Yn) + B(Yn, PnZ) + 
B(P n Z) 

(3.29) j t \\Y N f 2 + M 3 < c H^jv || 2 (IM 2 ~ 2A'i) + c \\Z\& , 

where I\\ is defined as Kq in ITTU1 but with a different c. We set 


cth 2 = inf jt £ (th 2 ,T) ||Yv(t)|| 2 > 2/vi j , 


Integrating (10)1) . we obtain 

||>Ar(o-H 2 )||2 + , 


\\Y n (t) H3 dt < \\Y N (m 2 


+ cM . 


Taking into account (TOTI) and choosing ||:ro|| 2 and M 2 sufficiently small, we obtain 

||^Ar(cr M2 ) Ha - 4 A Kl - 
It follows that cth 2 = T and that 
(3.30) ||^(T)|| 2 < 5, 

provided M and ||xo|| sufficiently small, which yields (1301) . 


4. Proof of Thfor.fm II .81 


As already explained, Theorem 11.81 follows from Proposition ll.151 We now prove 
Proposition HOI Let (a;J,a;g) £ (H 2 ) 2 . Let us recall that the process (A 1 , A 2 ) is 
defined at the beginning of section 3. 

Let <5 > 0, T £ (0,1) be as in Proposition 12.II and r defined in El, setting 


n = t, 


Tfc +1 = inf 


{ 


t > T k 


A 1 (^) 11 2 V ||A 2 (t)|| 2 < <5 j • 


it can be deduced from the strong Markov property of (A 1 , A 2 ) and from Propo¬ 
sition 130 that 

E(e" Tk+1 ) < K v E(e aTk (l + |A 1 ^)f + |A 2 (r fc )| 2 )) , 
which yields, by the Poincare inequality, 


{ E (e aTk+1 ) < cK” (1 + 2i5)E (e aTk ), 

E (e“ Tl ) < K” (l + \xl \ 2 + \x%\ 2 y 
It follows that there exists K > 0 such that 

E(e arfe )<A fc (l+|4| 2 + |x 2 | 2 ). 

Hence, applying Jensen inequality, we obtain that, for any 9 £ (0,1) 

(4.1) E ( e SaTk ) < K ek (l + |4| 2 + |ccg| 2 ) . 
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We deduce from Proposition 12.II and from (13.11) that 

¥(X 1 {T) ^ X 2 {T )) < 

provided are in the ball of (H 2 ) 2 with radius 5. 

Setting 

k 0 = inf {k £ N | X\r k +T)= X 2 (r k + T) } , 
it follows that fco < oo almost surely and that 

(4.2) P(fco>n)<Q) ■ 

Let 9 £ (0,1). We deduce from Schwartz inequality that 


E 


OO 

(e5“ T ‘ 0 ) = 


e f ar n 


Lfc 0 =n I < 


y/p (fc 0 > n) E(e e “ T "). 


Combining O and HD, we deduce 

E (ei“-o) < (jr (^0 j (l + \xl\ 2 + \x 2 \ 2 ) . 

Hence, choosing 9 £ (0,1) sufficiently small, we obtain that there exists 7 > 0 non 
depending on N £ N such that 

(4.3) E(e 7Tfc o)<4(l+|a;J| 2 + |^| 2 ). 

Recall that if (X^X 2 ) are coupled at time t £ TN , then they remain coupled for 
any time after. Hence X 1 ^) = X 2 (t) for t > r ko . It follows 

P (X^nT) y 2 X 2 {nT)) < Ae~ inT (l + |xj| 2 + |xg| 2 ) . 


Since (W 1 (nT), X 2 (nT)) is a coupling of {iV^ T )*8 x i,(V^ T )*5 x 2) : we deduce from 
Lemma m~n 


(4.4) 


(Kry S xh - (V&y f S xl 


< 4e 7 ” t ^1 + |xj|' + |x 2 | 2 ) , 


for any n £ N and any (xq,Xq) £ (H2) 2 . 

Recall that the existence of an invariant measure /tjv £ P(PnH) is justified in 
section 1.3. Let A £ P(H) and t £ R + . We set n = \Jp\ and C = 4e 7T . Integrating 
{xl,xl) over (( V t N _ nT )* A) ® hn in 14.41) . we obtain 


(V t N )*\- m <Ce-^(l+ [ \x 0 \ 2 \(dx)) 

var \ J H J 


which establishes 11.1911 . 


Appendix A. Proof of 12.201) 

For simplicity in the redaction, we omit 9 and N in our notations. 
Remark that 


(A.l) 


J = (VE (g(X(T))i/, x ) ,h) = J 1 + 2 J 2 , 
229 












Exponential mixing for the 3 D stochastic Navier-Stokes equations 


where 

l Ji = E((Vg(X(T)), V (T,0)-h)i/> x ), 

I J 2 = E [g{X(T)W x f 0 T ( AX(t ), A(r/(t, 0) • h)) ds) . 

According to 129, let us denote by D S F the Malliavin derivative of F at time s. We 
have the following formula of the Malliavin derivative of the solution of a stochastic 
differential equation 

D s X(t) = 1 t> a ri(t, s) ■ <f>(X(s)), 


which yields 


(A.2) 

where 


f D s X(t ) • m(s) ds = G(t) ■ to, 

J o 

G(t) ■ m = f rj(t, s) ■ </>(X(s)) ■ m(s) ds. 

J o 


The uniqueness of the solutions gives 


which yields 


Setting 

we infer from (IA.2I 


T}(t, 0) • h = r}(t, s) ■ ( r/(s , 0) • h) for any 0 < s < t, 

vi T , 0) • h = ^ J r}(t, s ) • (r?(s, 0) • h) ds. 
w(s) = cj)~ 1 (X(s)) ■ Tj(s,0) ■ h, 


(A.3) r,(T,0)-h=^G(T)-w=^J o D s X(T)-wds, 

which yields 

(Vg(X(T)), r/(T, 0) • h) = i J (V 5 (X(T)), D S X(T) ■ w) ds. 


Remark that 
It follows 


(■ D s g(X(T )), w) = (V. 9 (X(T)), D s X(T) ■ w ). 

(Xg(X(T)), g(T, 0) • h) = i J ( D s g(X(T )), w) ds, 

which yields 

(A.4) Ji = ^eJ %j) X (D s g(X(T)),w) ds. 

Recall that the Skohorod integral is the dual operator of the Malliavin derivative 
(See 2?j|). It follows 

(A-5) Jr = 1 e ^(A(T)) £ ^ x {w{i), dW^ . 

Recall the formula of integration of a product 

(A.6) f %l>x{w(t),dW(t)) = ipx [ (w(t),dW(t)) - f (D s i/j x ,w(s)) ds. 

Jo Jo Jo 
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Remark that 


D s ^ x = 2V4 [ AD s X{t) ■ ( AX(t))dt , 

Jo 


which yields, by {AX(t), AD s X(t) ■ w(s)) = (w(s), AD s X(t) ■ ( AX(t ))), 


rT r T 


/ (D s ip x ,w(s)) ds = 2tp x 

Jo 

We deduce from (TOli that 


(AX(t), AD s X(t) ■ w(s)) dtds. 


’ o jo 


r T 

(D s 'ijj x ,w(s)) ds = 2ip' x / t (AX(t), Ari(t,0) ■ h) dt. 
Jo 


(A.7) 

Remark that 

ipx = V’x = 0 if a < T. 
Hence combining (IA.6I) and CO- we obtain 

r T 


pCT p(T 

il>x(w(t),dW(t)) = %p x / {w(t),dW(t)) - 2 i/j x t {AX(t), Arj(t, 0) • h) dt. 
Jo Jo 


Thus, (12.20|) follows from (ED and 


References 

[1] J. Bricmont, A. Kupiainen and R. Lefevere, Exponential mixing for the 2D stochastic Navier- 
Stokes dynamics, Commun. Math. Phys. 230, No.l, 87-132, 2002. 

[2] P. Constantin and C. Foias, Navier-Stokes Equations, University of Chicago Press, Chicago, 
IL, 1988. 

[3] G. Da Prato, A. Debussche, Ergodicity for the 3D stochastic Navier-Stokes equations J. Math. 
Pures Appl. (9) 82 (2003), no. 8, 877-947. 

[4] G. Da Prato and J. Zabczyk, Stochastic equations in infinite dimensions, Encyclopedia of 
Mathematics and its Applications, Cambridge University Press, 1992. 

[5] G. Da Prato and J. Zabczyk, Ergodicity for Infinite Dimensional Systems, London Mathe¬ 
matical Society Lecture Notes, n.229, Cambridge University Press, 1996. 

[6] A. Debussche, C. Odasso, Ergodicity for the weakly damped stochastic Non-linear Schrodinger 
equations, Journal of Evolution Equations , Birkhuser Basel, vol. 5, no. 3, pp. 317-356. 

[7] A. Debussche, C. Odasso, Markov solutions for the 3D stochastic Navier Stokes equations 
with state dependent noise, Journal of Evolution Equations, Birkhuser Basel, vol. 6, no. 2, pp. 
305-324. 

[8] W. E, J.C. Mattingly, Y. G. Sinai, Gibbsian dynamics and ergodicity for the stochastically 
forced Navier-Stokes equation, Commun. Math. Phys. 224, 83-106, 2001. 

[9] F. Flandoli Irreducibilty of the 3-D stochastic Navier-Stokes equation, Journal of Functional 
Ananiy sis, 149, 160-177, 1997. 

[10] F. Flandoli, An introduction to 3D stochastic Fluid Dynamics, CIME lecture notes 2005. 

[11] F. Flandoli and D. Gatarek, Martingale and stationary solutions for stochastic Navier-Stokes 
equations, PTRF, 102, 367-391, 1995. 

[12] F. Flandoli and M. Romito, Partial regularity for the stochastic Navier-Stokes equations, 
Trans. Am. Math. Soc. 354, No 6, 2207-2241, 2002. 

[13] F. Flandoli and M. Romito, Markov selections for the 3D stochastic Navier- Stokes equations, 
in preparation. 

[14] Foias, C. and Temam, R., Gevrey class regularity for the solutions of the Navier-Stokes 
equations, J. Funct. Anal. 87 (1989), no. 2, 359-369. 

[15] M. Hairer, Exponential Mixing Properties of Stochastic PDEs Through Asymptotic Coupling, 
Proba. Theory Related Fields, 124, 3 :345-380, 2002. 

[16] M. Hairer, J. C. Mattingly Ergodic properties of highly degenerate 2D stochastic Navier- 
Stokes equations C. R. Math. Acad. Sci. Paris 339 (2004), no. 12, 879-882. 

[17] M. Hairer, J. Mattingly, Ergodicity of the 2D Navier-Stokes equations with degenerate forcing, 
preprint. 


231 



Exponential mixing for the 3 D stochastic Navier-Stokes equations 


[18] W.D. Henshaw, H.O. Kreiss, L.G. Reyna, Smallest scale estimate for the Navier-Stokes 
equations for incompressible fluids, Arch. Rationnal Mech. Anal. 112 (1990), no. 1, 21-44. 

[19] G. Huber, P. Alstrom, Universal Decay of vortex density in two dimensions , Physica A 195 , 
448-456, 1993. 

[20] Karatzas I. and Shreve S. E., Brownian motion and stochastic calculus, Second edition, 
Graduate Texts in Mathematics, 113. Springer-Verlag, New York, 1991. 

[21] S. Kuksin, On exponential convergence to a stationary mesure for nonlinear PDEs , The M. 
I. Viishik Moscow PDE seminar, Amer. Math. Soc. Trans. (2), vol 206, Amer. Math. Soc., 
2002 . 

[22] S. Kuksin, A. Shirikyan, Stochastic dissipative PDE’s and Gibbs measures , Commun. Math. 
Phys. 213 , 291-330, 2000. 

[23] S. Kuksin and A. Shirikyan, Ergodicity for the randomly forced 2D Navier-Stokes equations, 
Math. Phys. Anal. Geom. 4, 2001. 

[24] S. Kuksin, A. Shirikyan, A coupling approach to randomly forced randomly forced PDE’s I, 
Commun. Math. Phys. 221, 351-366, 2001. 

[25] S. Kuksin, A. Piatnitski, A. Shirikyan, A coupling approach to randomly forced randomly 
forced PDE’s II, Commun. Math. Phys. 230 , No.l, 81-85, 2002. 

[26] S. Kuksin, A. Shirikyan, Coupling approach to white-forced nonlinear PDEs, J. Math. Pures 
Appl. 1 (2002) pp. 567-602. 

[27] J. Mattingly, Exponential convergence for the stochastically forced Navier-Stokes equations 
and other partially dissipative dynamics, Commun. Math. Phys. 230 , 421-462, 2002. 

[28] J. Mattingly, E. Pardoux, Ergodicity of the 2D Navier-Stokes Equations with Degenerate 
Stochastic Forcing , preprint 2004. 

[29] Nualart, Malliavin Calculus and related topic, Probability and its Applications, 1995, 
Springer. 

[30] C. Odasso, Ergodicity for the stochastic Complex Ginzburg-Landau equations, to appear in 
Annales de l’institut Henri-Poincare, Probabilites et Statistiques. 

[31] C. Odasso, Exponential Mixing for Stochastic PDEs: The Non-Additive Case, preprint avail¬ 
able on http://www.bretagne.ens-cachan.fr/math/people/cyril.odasso/ 

[32] A. Shirikyan, Exponential mixing for 2D Navier-Stokes equation pertubed by an unbounded 
noise, J. Math. Fluid Mech. 6, no. 2, 169-193, 2004. 

[33] A. Shirikyan, Analyticity of solutions of randomly perturbed two-dimensional Navier-Stokes 
equations, Uspekhi Mat. Nauk 57 (2002), no. 4(346), 151-166; translation in Russian Math. 
Surveys 57 (2002), no. 4, 785-799. 

[34] R. Temam, Navier-Stokes Equations. Theory and Numerical Analysis., North-Holland, 
Amsterdam-New York-Oxford,1977. 

[35] R. Temam Navier-Stokes equations and nonlinear functional analysis, Philadelphia (PA US) 

: SIAM , 1995 , CBMS-NSF regional conference series in applied mathematics ; 66 ISBN 
0-89871-340-4 

[36] R. Temam, Some developments on Navier-Stokes Equations in the second half of the 20th 
Century, ”in Developpement des Mathematiques au cours de la seconde moitie du XXeme 
siecle”, J.P. Pier ed., Birkhauser Verlag, Basel, 2000. 


232 



